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!— I ■ Abstract 

In this paper, we investigate the dynamical formation and evolution of 2 + 1-dimensional 



(N 
> 



charged black holes. We numerically study dynamical collapses of charged matter fields in 
an anti de Sitter background and note the formation of black holes using the double-null 
' formalism. Moreover, we include re-normalized energy-momentum tensors assuming the S- 

wave approximation to determine thermodynamical back-reactions to the internal structures. 
If there is no semi-classical effects, the amount of charge determines the causal structures. If 
the charge is sufficiently small, the causal structure has a space-like singularity. However, as 
the charge increases, an inner Cauchy horizon appears. If we have sufficient charge, we see a 
space-like outer horizon and a time- like inner horizon, and if we give excessive charge, black hole 
i horizons disappear. We have some circumstantial evidences that weak cosmic censorship is still 

satisfied, even for such excessive charge cases. Also, we confirm that there is mass inflation along 
the inner horizon, although the properties are quite different from those of four-dimensional 
cases. Semi-classical back-reactions will not affect the outer horizon, but they will affect the 
inner horizon. Near the center, there is a place where negative energy is concentrated. Thus, 
charged black holes in three dimensions have two types of curvature singularities in general: via 
mass inflation and via a concentration of negative energy. Finally, we classify possible causal 
structures. 
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1 Introduction 



Traditionally, people regarded that three-dimensional models are good toy models with which to 
study the basic nature of pure gravity and quantum gravity [1]. First, there may be a qualitative 
relation between a cosmic string in four dimensions and a point source in three dimensions [2J. 
Second, three-dimensional gravity is simpler than the four-dimensional version, making it useful 
for obtaining some intuitions about quantum gravity using gauge theoretical correspondences [3]. 
However, given that this is too simple, it was doubtful that a black hole solution would be found 
in a three-dimensional background. 

The first analytic solution in three-dimensional gravity was found by Banados, Teitelboim, and 
Zaneli [4] [5]. They found the black hole solution in an anti de Sitter background. In their solutions, 
the negative cosmological constant played a crucial role. These solutions can have mass, charge, 
and angular momentum. Moreover, the solutions can have the Hawking temperature [6|. However, 
they are thermodynamically stable because the specific heat is positive [7]. 

The so-called BTZ solutions are useful in the context of holography. According to the AdS/CFT 
correspondence [5], there is a relationship between the bulk three-dimensional gravity and the 
conformal field theory at the asymptotic two-dimensional boundary of such a three-dimensional 
bulk. That is, we may relate the gravitational perturbation on the bulk to the thermal excitation 
on the boundary field theory. Moreover, it is known that the U-duality transformation maps four 
or five-dimensional black holes to three-dimensional BTZ black holes [9]; here, these black holes, 
which are connected by the U-duality, have the same number of microstates. Hence, the microscopic 
derivation of the entropy can be replaced by the study on the corresponding derivation for the BTZ 
black hole [10] . In these senses, three-dimensional black holes are more important than just an 
artificial toy model. 

In this paper, we will begin to study the dynamics of 2 + 1-dimcnsional charged black holes, 
including gravitational collapses and back-reactions, using the semi-classical effects of Hawking 
radiation. We solved time-dependent Einstein and field equations using the double-null formalism 
[TT j [12] P ~5 ] ■ Of course, it is too ambitious to extend this method to holography or string theory 
in this paper. However, if we establish a method with which to study dynamical black holes 
in an anti de Sitter background, it will eventually become possible to study the time-dependent 
correspondences between the bulk and boundary someday. 

Thus far, there is less reference to compare our calculations on the dynamic formation of three- 
dimensional charged black holes (e.g., [2]), whereas there is some literature on the gravitational 
collapses of neutral matter [15] , especially in the context of critical collapses [16] . On the other hand, 



3 



in four dimensions, there are numerous studies of gravitational collapses and evaporations |llj | 12 ) . 
One important feature in four-dimensional dynamical charged black holes is the existence of mass 
inflation |17j . Mass inflation is a phenomenon in which the local mass function of an observer near 
the Cauchy horizon increases exponentially. If it is true, then the internal structure of charged 
black holes cannot be described by the Reissner-Nordstrom solution |18] and the internal structure 
should be significantly modified [19] . This was directly confirmed by numerical simulations, and 
there are more related works |ll)|12j. Three-dimensional gravity will be similar as regards this 
point, as it is expected that there is mass inflation in general in charged black holes 20 . However, 
there are some expectations that three-dimensional mass inflation is different from four-dimensional 
cases [21] , Moreover, because black holes are thermodynamically stable in three-dimensions, it will 
be interesting to determine what the semi-classical back-reactions are in the internal structures of 
black holes (there are discussions on the semi-classical effects in the static limit [22 ). 

This paper is organized as follows: in Section [2j we summarize analytic solutions for three- 
dimensional charged black holes and discuss the mass inflation issue. The paper will then provide 
the theoretical motivations to study three dimensions. In Section [21 we describe our model of the 
numerical simulations, the implementation to the double-null formalism, and comment on the initial 
conditions and dimensional analysis. In Section^ we report the results. First, we report results for 
cases without semi-classical effects. Second, we include semi-classical effects and observe how they 
modify causal structures. Finally, in Section [5j we conclude and comment on causal structures and 
cosmic censorship issues. 

2 Preliminaries 

In this section, first we briefly review the analytic and static charged black hole solutions of three 
dimensions. Second, we summarize thermodynamic properties and internal stability issues via mass 
inflation regarding their solutions. 

2.1 Static solutions: causal structures and thermodynamics 
2.1.1 Three dimensions: BTZ solutions 

For three dimensions, the first simple black hole solution in the anti de Sitter background is known 
to be the BTZ solution (for non- rotating cases) [4]. For the Einstein-Maxwell action 




(1) 
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the solution becomes 

ds| TZ = -N 2 dt 2 + N~ 2 dr 2 + r 2 dip 2 , (2) 

where 

N 2 (r)=~M+ r ^-Q 2 \n r ^ : (3) 

A = — l/l 2 , Q is electric charge, and M is interpreted by the mass. 

If Q = 0, then we obtain a neutral black hole solution with r + = lyM. Thus although this is 
in the M = limit, the center still contains a singularity (Q = M = in Figure [1]). The proper 
anti de Sitter limit is obtained by the M = — 1 limit: 



ds 2 = - ( 1 + — )df 2 + 1 + — ) dr 2 + r 2 dip 2 . (4) 



If Q =/= 0, then there can be two horizons: the outer horizon r + and the inner horizon r_. If 
M > 1, then there are always two horizons with r + > iVM and r_ < l\f~M (Case 1 in Figure [lj. 
However, if M < 1, then the number of horizons depends on Q: for a small Q limit, there can be 
two horizons lyM > r + > r_ (Case 1 in Figured]); for a large Q limit, there is no horizon (Case 3 
in Figure [lj; and the extreme limit, r± = Ql and M — Q 2 (l — InQ 2 ) (Case 2 in Figure [TJ. In 
the neutral limit, the singularity r = is space-like (Q = in Figure [TJ; in the charged limit, the 
singularity is time-like (except Q = M = 0). Therefore, it is not difficult to draw Penrose diagrams 
in the static limit [S]. 

The Hawking temperature Tbtz becomes 



1 dN 2 

1 BTZ 



47r dr 



" + ^-^], (5) 



and Tbtz is positive or zero because r + is greater than the extreme limit size Ql [I] [B] . The entropy 
Sbtz becomes 

Abtz irr + 

i>BTZ - = (o) 

However, the specific heat Cbtz becomes 

->2;2 



SSbTZ (, , Q 2 l 



and hence the black hole is in thermal equilibrium and the outer horizon does not shrink [7]. 
Therefore, when we consider the dynamical formation of a black hole, the outer horizon should be 
always space-like or null. 

Note that in the double-null coordinates, 

c?s B tz — —a 2 dudv + r 2 d(p 2 , (8) 
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Case 1 Case 3 

Figure 1: Penrose diagrams of BTZ solutions for Q = and M > 0, M = Q = 0, Case 1 (there 
are two horizons), Case 2 (there is one horizon), and Case 3 (there is no horizon but a time-like 
singularity) . 
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the Einstein equation for the vv component around the outer apparent horizon r t = is 

v ,w — &*kT vv . (9) 

Therefore, the outer apparent horizon is space-like or null if and only if T vv > 0. In other words, 
for BTZ solutions, even when they are dynamically generated with a Hawking temperature, T vv 
should be positive or around the outer apparent horizon. 

2.1.2 Instability of Cauchy horizons 

There are a number of discussions on the instability of Cauchy horizons r_ [3D] [21] [23] • We can 
re-write the metric in the following form: 

ds 2 = -N 2 dv 2 + 2dvdr + r 2 dQ 2 D _ 1 . (10) 

Here, D is the space dimensions (2 or 3). Now, without loss of generality, we can ignore the angular 
part to discuss in-going or out-going null geodesies. 

Thus, for an in-falling matter along the in-going null direction (for coordinates [v, r]), the energy- 
momentum tensor components T a p can be represented by 

F(v) 

Tap oc -p—fd a vdpv, (11) 

where F(v) is an arbitrary function that refers to the luminosity function |20j |21j. At this point, 
if there is an observer who moves along the out-going null direction and approaches the Cauchy 
horizon r_, then for a null geodesic of the observer P, (d a v)l a — l v — dv/drj and dv/drj oc exp/i_u 
around the v — > oo limit, where 77 is an afhne parameter of the observer and k_ is the surface 
gravity of the inner horizon. In conclusion, the out-going observer who approaches the Cauchy 
horizon measures the energy density for the energy flow that flows along the Cauchy horizon by 

mm- 

p = T ap l a f oc ^e 2K - v . (12) 

Note that for four dimensions, for a realistic black hole the behavior of the luminosity function 
should be polynomial to F(v) ~ v~ p [24 . Although there is no such guidance regarding the 
luminosity function for three dimensions, if we believe that such a polynomial fall-off is reasonable 
for three-dimensional cases, then we can conclude that for both three- and four-dimensional cases, 
any small energy flow along the Cauchy horizon is amplified for an observer who approaches to the 
Cauchy horizon. This effect is known as mass inflation. 

Therefore, the inner horizon cannot be stable and there should be strong back-reactions due to 
high energy density [19) . In many cases, it is believed that the Cauchy horizon should be a curvature 
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singularity |25j . Hence, we cannot trust the static solutions to study the inside structures of charged 
black holes. Therefore, it is necessary to introduce numerical tools to study these structures. 

2.2 Motivation of this paper 

We have two main motivations for this paper. 

1. Theoretical interests: Three-dimension is important. First, it is a toy model of a black 
string that may be possible in the early universe. Second, as we commented in the introduc- 
tion, in string theory, there is a holographic correspondence between bulk black holes in an 
anti de Sitter space and its boundary conformal field theory. One can relate the gravitational 
perturbation on the bulk to the thermal excitation on the boundary field theory and one 
may restrict the perturbation by the spherically symmetric scalar sector. Although the effect 
of the scalar perturbation on the boundary field theory is not yet revealed fully, our study 
and methodology on dynamical black holes may shed light on the study of the dynamical 
correspondence between gravity and field theory. 

2. Numerical interests: According to previous discussions, charged black holes in three dimen- 
sions will have mass inflation around the inner horizon. This means that the inner horizon and 
internal structures cannot be described by analytical methods. To induce this mass inflation, 
the black hole should be fully dynamical. Then it is inevitable to use numerical methods. 

Therefore, there are various interests to study three-dimensional black holes. To study internal 
structures, it is inevitable to use the numerical methods. In this paper, we study the dynamical 
formation of spherical symmetric charged black holes and its back-reactions for three dimensions. 
We prepare a complex scalar field with the Maxwell field; it is used to form a charged black hole. 
Then, it will make it possible to understand mass inflation and the possible curvature singularities 
inside of three-dimensional black holes. 

Moreover, we can include Hawking radiation by introducing re-normalized energy-momentum 
tensors with the S'-wave approximation. We can check whether the horizon behavior and the energy 
condition are consistent with the thermodynamics results, even when we include the re-normalized 
energy-momentum tensors. We can also determine the back-reaction of the re-normalized energy- 
momentum tensor to the inside structure of the black hole. 
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3 Model for three-dimensional gravity 

In this section, we introduce a model for three-dimensional charged black holes. To use numerical 
simulations, we introduce the double-null formalism and discuss the details of the initial conditions 
and free parameters for the simulations. Also, we briefly comment on the dimensional analysis of 
each functions and parameters. 

3.1 Theory 



Let us introduce the action in three dimensions by 



/ 



S = I dx 3 ^ 



1 ' (R - 2A) - \ + ieA^) gT (</>. v - ieA v ~$) - r^-F lu ,F» u 



(13) 



_16tt v ' 2 " T ' 16tt 

where R is the Ricci scalar, is a complex scalar field with a gauge coupling e and a gauge field 
A^, where F^ = A v;ii - A^ v . 

The Einstein equation becomes as follows: 

Gfj, v + Ag^ = 87rT MI/ , (14) 

where the total energy-momentum tensors are 

V = TC+(f«). (15) 

The classical part Tf? v is 

T^V = \{^l ) -,»+~4>^ ( t>-M) 

+ \ {-<P;n ieA i<<P + (p-^ieA^ + <t>.^ieA v <j> ~ ^ieA^cf) 

+ ^-F^FS + e 2 A„i4„# + Cg^ (16) 

and (T^) is the re-normalized energy-momentum tensor part to include semi-classical effects, where 
C = -(1/2) + ieA^) g»» - ieA v J) - (1/16^)F^F^. 

The field equations for the complex scalar field </> and the Maxwell field A^ are as follows: 

+ ieA»(2<j>.f M + ieAf l <l>)+ieAf l . v g» v <l> = 0, (17) 
^F\, v - ie<t> - ieAfJ) + ie~$ + leA^) = 0. (18) 

3.2 Implementation via the double-null formalism 

We use the double-null coordinates 

ds 2 = —a 2 (u, v)dudv + r 2 (u, v)dif 2 , (19) 
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assuming spherical symmetry, where u is the retarded time, v is the advanced time, and ip is the 
angular coordinate. 

We follow the notation of [26] jTT] [T2] [13] : the metric function a, the radial function r, and a 
scalar field s = y/incj), and define 

h=—, d=—, f = r u , g = r yVl w = s, u , z = s v . (20) 
a a 

The Einstein tensors are then given as follows: 



G uu = --(/,„ -2//0, (21) 
r 

Guv = i (22) 
r 

G vv = —(g >v -2gd), (23) 
r 

4r 2 

G w = T d u . (24) 

Via the gauge symmetry of the Maxwell field, we can choose the Maxwell field — (a, 0, 0) by 
introducing the function a(u, v). Thus, we can obtain the energy-momentum tensors of the classical 
part with the cosmological constant term by: 

T uu ~ tt9uu = t~ [ww + iea(ws - ws) + e 2 a 2 ss] , (25) 

07T 47T 

T ^~^9u V = ^J + Y^A, (26) 

T vv - g^9vv = ~^. ZZ i ( 27 ) 



c A r 2 



(wz + zw) + tea(zs — zs) H - — 



8^ < 28 > 



^ 87r yw 47ra 2 

Note that F M „ :M = A-kJ v , where J v is asymptotically interpreted by the electric current. If 
we define q(u,v) = ra tV /a 2 as the charge function, then J„ = (l/27rr)(— q tU , q tV , 0) and hence 
corresponds to the concept of the charge. 

There are several discussions pertaining to re-normalized energy-momentum tensors for three 
dimensions in the literature [22]. As expected for two- or four-dimensional cases [23] [20], the re- 
normalized energy-momentum tensors diverge along the inner horizon for the static metric. How- 
ever, we know that the internal structures cannot be described by the static metric for fully dynam- 
ical cases. Due to dynamical situations and spherical symmetry, therefore, the following represents 
the best choice for the re-normalized energy-momentum tensor: use the (1 + l)-dimcnsional results 
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divided by 2irr and apply the S'-wave approximation by 



(t H J = ^(^~ /l2 )' ( 29 ) 

= = ~d, u , (30) 
Zirr 

(f*) = £;(dv-d 2 ), (31) 



where P = AfK/ 12tt and A/" is the number of massless scalar fields. 
After solving these coupled equations, we can write all equations: 

/,„ = 2A-8 7 rr(T u c „ + (fH)) , (32) 

g, v = 2gd-8nr(T^ + (f^, (33) 

9,u = f,v = ^r(TC-A guv ^ +8 KP^(T^-^g w y (34) 

A -h - -^L (tc _ A ^ ^ 

",u — "-,[1 — 2 I <P<P o_5w ] > l^ / 



8tt' 



including the field equations 



a 2 g 
r 



(36) 



q.v = -—(sz-sz), (37) 

/z gw . ieags ie 2 

z„ = «) t = zeetz a as. (38 

2r 2r 2r 2r ^ V ' 

At this point, the equations of a tUV , r utl , and s jM „ parts can be represented by first-order 

differential equations. We can then implement the same integration scheme used in previous papers 

[T2] [IH]- We used the second-order Runge-Kutta method [30]; however, due to the differentials 

in Equations (l32|) and (j33l) . our results will converge to the first order if P > 0. Tests of the 

convergence are provided in Appendix B. 

3.3 Definition of mass function 

One important question when interpreting our results is whether or not these black hole solutions 
are BTZ solutions. Of course, due to the Birkhoff theorem for three-dimensional cases, we can be 
sure that our solutions will be BTZ (for further discussions, see Appendix A). Hence, there will be 
a natural concept for the mass function. 

First, we can define the general spherical symmetric metric in {D + 1) dimensions in the static 
limit by the following coordinates: 

ds 2 = + +r > d£l 3 D _ 1 (39) 

= -N 2 dUdV + r 2 {U,V)dti 2 D _ 1 , (40) 



11 



where 



dU = dt 
dV = dt 



dr 

dr 
7V2- 



(41) 
(42) 



In many cases, function N vanishes near the horizon; hence, it will be convenient to redefine the 
coordinate in the following form: 



ds 2 = —a 2 (u,v)dudv + r 2 (u,v)dti 2 D _ ll 
where the coordinate transformation is 



dr 
dt 



Thus, we can show 



N 



r^ u du + r iV dv, 
a 2 / dv du\ 
4 V r ,u r ;V ) 



(43) 

(44) 
(45) 

(46) 



Note that r tU r tV is negative outside of the black hole. Therefore, the signature is fine. Moreover, 
we can easily see that the apparent horizons r v = or r „ = make N = 0, although a is not zero. 
Thus, in the four dimensions, we can define the mass function |31] by 



r I . q 

m {i) {u,v) = - + + - 



2 A , 



Also, for three dimensions, we can similarly define the mass aspect by 



mr 3) (u,v) 



q 2 In- 



4?* vT ,v 



(47) 



(48) 



This definition is consistent with the findings in the literature |16j . 
3.4 Initial conditions and free parameters 

We need initial conditions for all functions (a,h,d,r,f,g,s,w,z,a,q) on the initial u — u\ and 
v = Vi surfaces, where we set ui = V\ = (Figure [2]). 





Fixed by hand 


Fixed by constraints 


Fixed by evolution equations 


V = Vi 


s = 0,/,r 


w = h = a = q = 0, a = 1 


d,g,z 


U = Ui 


s,g,r 


z = s tV , d = 2rzz, a <s— d = a, v ja 


h,f,w,q,a 



Table 1: Summary of the assignments of initial conditions for v = v\ (the in- going initial surface) 
and u = Ui (the out-going initial surface) . 
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Figure 2: Integration domain of our simulations. We assign initial conditions along the in-going 
and out-going null surfaces along the arrows. 

We have gauge freedom to choose the initial r function. Although all constant u and v lines 
arc null, there remains freedom to choose the distances between these null lines. Here, we choose 
r(0, 0) = 7*0, f(u, 0) = r u o, and g(0,v) = r v o, where r U Q < and r„o > such that the radial 
function for an in-going observer decreases and that for an out-going observer increases. 

We use a shell-shaped scalar field. Therefore, its interior is not affected by the shell. First, it 
is convenient to choose r u o = —1/2 and r v o = 1/2; the mass function on u\ = v\ = then becomes 
m( 3 )(iti, «i) = Tq/1 2 — l/a(u- u t>i) 2 , as initially there is no charge. The pure anti de Sitter limit is 
r7i(3) — — 1) an d hence as long as 777( 3 ) > —1, the choice of a(ui,Vi) is physically allowed, as one 
can always assume that there is a certain core of mass at the center. In this paper, we choose 
a(u, 0) = 1 for convenience. Also, s(u, 0) = w(u, 0) = h(u, 0) = a(u, 0) = q(u, 0) = hold. 

We need more information to determine d, g, and z on the v = surface. We obtain d from 
Equation fl35]), g from Equation (|3"3)). and z from Equation fl^B} (Tabled]). 

We can choose an arbitrary function for s(0, v) to induce a collapsing pulse. In this paper, we 

use 

s(ui,v) = A sin 2 ( 7r ) exp ( 27T7 ) (49) 

\ Vf-ViJ V Vt-ViJ 

for Vi < v < Vf and otherwise s(u{, v) = 0, where Uf = 0, Vi = 0, and Vf = 20 denotes the end of the 
pulse in the initial surface. We then obtain z(0, v) = s(0, v) jV . This implements one pulse of energy 
(T vv ~ z 2 ) along the out-going null direction by the continuous function z(0,v). 

Also, from Equation (1331) . we can use d = Irzz on the u = surface, given the assumption that 
there is no Hawking effect on the initial surface. We then obtain d{Q,v). By integrating d along v, 
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we have a(0, v). 

We need more information for h, /, q, a, and w on the u — surface. We obtain h from Equa- 
tion (1351) . / from Equation ([34|), a from Equation p6[) . g from Equation p7p . and it) from Equa- 
tion P"5)l . This finishes the assignments of the initial conditions (Table [TJ. 

We choose ro = 10, A/87T = —0.0001, and A = 0.2, leaving the two free parameters (e, P), where 
e is the gauge coupling and P determines the strength of the semi-classical effects. 

3.5 Dimensional analysis 

For four-dimensional cases, if we define c = G4 = 1, then the Planck units (Planck length Zpj^ and 
Planck mass mpi^) become 

IpiA = TlP1.4 = Vh. (50) 

Because P oc Afh should be fixed by a simulation parameter, as Af becomes larger and larger, h 
becomes smaller and smaller. For example, for four-dimensional simulations, if we made a black 
hole with size r + = 10, we fixed P = 1, and we assume Af = 100, then the physical size of the 
black hole in Planck units is r + = 100 ipi,4. By the same argument, the 'bare' simulation results 
with dimension length L, mass M, or curvature R (e.g., the Ricci scalar) are re-scaled in the Planck 
units by 

L = VAfLl P i A , (51) 
M = VNMmpiA, (52) 
R = ^ P1 2 4 , (53) 

where Afh = Afl% lA = Afmp 1A = 1 [32]. 

For three-dimensional cases, the dimensional analysis is different. If c = G3 = 1, then 

lpi,3 = fi (54) 

and mpi^ is dimensionless. Again, P can be fixed as a constant (e.g., 1), and then any bare 
quantities with dimensions length L, mass M, or curvature R are re-scaled in Planck units by 

L = 7ViZpi )3) (55) 
M = Mmpi, 3 , (56) 

where Afh = Aflpi 3 = 1. 
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One interesting difference between the three-dimensional and four-dimensional results is the 
large- TV cutoff. According to Dvali [33], when there is a large number TV of particle species that 
contribute in the mode of Hawking radiation, the cutoff scale should be stretched by 

^cutoff,4 = VA/7pi,4- (58) 

Dvali supported this idea using arguments that included black hole complementarity |34j : if a 
black hole is smaller than / cu toff,4, then the evaporation is too fast and one cannot regard it as a 
semi-classical black hole. In such a case, it can violate the principle of black hole complementarity^ 

In the double-null formalism, for four-dimensional cases the Dvali cutoff is naturally implemented 
via the (1 - P/r 2 )- 1 term in an equation |12j . Thus, \[P ~ V^V7pi,4 was the natural cutoff scale 
of the simulation. Therefore, any results of the double-null simulation was consistent with Dvali's 
argument. However, in the three-dimensional case, there is no singular term despite the fact that 
we can choose a large P parameter. It is well known that a black hole is in the thermal equilibrium 
and does not totally evaporate in three-dimensions [7]; hence, we do not have to worry about the 
information loss problem. Therefore, there is no essential reason to consider a large-TV cutoff. This 
thermodynamic intuition is naturally implemented in our double-null formalism. 

On the other hand, the correction term from the re-normalized energy-momentum tensor for 
four dimensions becomes 

(f^ v ) - ^(dimensionless terms) + 0{h 2 ) + (59) 

while for three-dimensional cases, 

(TH) ~ — (dimensionless terms) + 0(h 2 ) + ... (60) 

occurs. Therefore, r = \fP in four dimensions and r = P for three dimensions are a type of 
turning point during which the contribution of the re- normalized energy-momentum tensor becomes 
dominant over the other terms. If this intuition is correct, the inside of the r = P region for three- 
dimensional cases is dominated by negative energy via semi-classical effects rather than classical 
matters. This may imply that r < P is a physically special region that changes the classical 
properties, although it is not a definite cutoff as in four-dimensional cases. 

1 Also, see further discussions on black hole complementarity in |35| . 
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4 Dynamical formation and evolution of three-dimensional 
black holes 



4.1 Gravitational collapses without Hawking radiation 

We study the gravitational collapse of a charged matter shell and the formation of black holes 
without Hawking radiation. We varied the gauge coupling parameter e to change the amount of 
charge of the matter shell. Figure [3] is the result of the numerical simulations. 

4.1.1 Causal structures 

For the e = limit, there is no gauge coupling and the black hole is neutral. There is a space- like 
singularity and an apparent horizon. The apparent horizon grows in the space-like direction during 
the matter collapse (v < Vf = 20) and approaches a null direction (Type 1 in Figure [S]). 

Also, if gauge coupling is sufficiently small so that the amount of charge is negligible (e = 0.075), 
the causal structure is similar to the neutral case. However, the location of the singularity changes 
slightly. As we increase the parameter e, the location of the space-like singularity drastically changes 



If we increase e more, we cannot see the space-like singularity (e = 0.2). The result of e = 0.2 
itself does not offer information about the existence or absence of the space- like singularity. However, 
in the e = 0.3 and e = 0.38 cases, we can see a time-like piece of the apparent horizon during the 
matter collapse. Although the apparent horizon is time-like, it does not violate the null energy 
condition, as an out-going observer sees a decreasing r. Rather, the time-like part should be 
identified with a type of inner horizon of the black hole; the outer and inner horizons are separated 
during gravitational collapses. Thus, although we cannot see beyond the u — 20 surface, we can 
guess that the inner horizon should be time-like. Also, we believe that there should be a singularity 
around the center via the singularity theorem, and that it may be inside of the inner horizon. If 
this expectation is reasonable, then the singularity should be time-like. 

Now let us comment on the behavior of r „ with the large v limit. In Figure El we plot the 
behaviors of log|r it ,| as a function of \ogv along u = 10. Note that for e = 0.15, 0.2, and 0.38, 
there are horizons and the u — 10 surface will be included by the horizon, while the u = 10 line for 
e = 0.4 is not trapped by a horizon. The plot shows that |r „| converges for a large v limit as a 
function of , where /3 > 1 is approximately 2.3, 2.5, 1.6, and 3.8 respectively for the four cases. 
Note that if /3 is larger than 1, then 



(e = 0.15). 




(61) 
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Figure 3: Numerical results for P = (without Hawking radiation) while e varies. 
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Figure 4: log|r„| versus logw along it = 10 for e = 0.15, 0.2, 0.38, and 0.4. The gradients for a 
large v limit are -2.2951, -2.4932, -1.6447, and -3.7969, respectively. 

will converge for v — > oo limit. This implies that, in the v — > oo limit, the inside of the horizon 
becomes a null r „ = horizon with a finite r. This behavior arises even for the e = 0.4 case, as it 
is in the anti de Sitter background. In this case, a large v limit should be regarded as the time-like 
boundary of the anti de Sitter space, implying that an outgoing observer cannot experience an 
increase of the area. 

Therefore, in terms of the inside of the outer apparent horizon, as in the four-dimensional cases, 
given the large v limit, we can interpret these as an inner Cauchy horizon (r„ — > 0). In this case, 
the inner horizon is null. Therefore, we can determine the causal structure for a small charge case 
(Type 2 in Figure O and for a large charge case (Type 3 in Figure [5]). 

If we increase e until it is excessively large, we cannot see apparent horizons and the black hole 
disappears. This phenomenon appears odd. For four-dimensional cases, we can make an asymptotic 
matter shell with M < Q, but a black hole always forms for any field combinations [12]. This means 
that, for four-dimensional cases, the excessive charge is scattered and repulsed via electromagnetic 
forces (by purely classical effects). This makes it impossible to see a naked singularity (Q > 
M of Reissner-Nordstrom solutions) and conserves weak cosmic censorship. However, for three- 
dimensional cases, we cannot not see a horizon in this case. To interpret this, there are two 
possibilities: a naked singularity (Type 4-1 in Figure [3]) or a regular center (Type 4-2 in Figure [5]). 
In principle, it is impossible to decide for certain, but we can offer a reasonable guess as to what 
the correct answer is, as shown in Section [4.1.21 



IS 



singularity 




Type 4-1 Type 4-2 



Figure 5: Causal structures without Hawking radiation. Type 1: Q = case. Type 2: if the 
charge is sufficiently small, there can be a space-like singularity as well as a null inner Cauchy 
horizon in the v oo limit. Type 3: if the charge is sufficiently large, one cannot see a space-like 
singularity. The singularity is then beyond the inner horizon and it must be time-like; therefore, it 
is a reasonable to guess that there is a time-like inner horizon beyond the Cauchy horizon. Type 4-1 
and Type 4-2: if the charge is excessive, the horizon then disappears. In this case, there may be a 
time-like singularity beyond the Cauchy horizon. Or, the charged matter may be repulsed via the 
strong electric field and fail to form a singularity. 
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4.1.2 Black hole formation near extreme limit 

We observe the details of causal structures between e = 0.35 and e = 0.4, where a type of critical 
behavior occurs (Figure 15]). As we increase the charge, two horizons (time-like inner and space- like 
outer apparent horizons) approach. Additionally, the advanced time v when the apparent horizon 
appears seems to increase so that the large charge postpones the formation of the apparent horizon. 
The location of the apparent horizon decreases, which implies that the effective mass of a black 
hole decreases as the charge increases. 

We know that there will not be horizons when e = 0.4 or e = 0.5 between v = and v = 180, 
even beyond the u — 20 line (we know that the horizons disappeared as the charge increases, as 
shown in Figure [SJ . Therefore, if there is a singularity for such cases, it should imply that it is a 
type of naked singularity. If this is true, then it is fair to say that we cannot see beyond the u = 20 
line in principle. However, reasonable assumptions can be made given some circumstantial evidence 
and extrapolations. 

We plot charge q as a function of r at the advanced time v — 180 in Figure [7J As the charge 
increases, the charge of outside (large r) increases while the charge of inside (small r) decreases. 
This implies that our field combination asymptotically induces a large charge, as we use a large 
coupling e. However, as the charge increases, mutual repulsion is noted. Hence, a collapse near the 
center is not probable, and the center has less charge as e increases. 

Note that this diagram is a log-log plot and that appears to be linear for a large r and a small 
r. Of course, we are interested in the charge distribution near the center, and hence, an estimation 
of the gradient of the logg-logr plot near the center, i.e., the exponent of q = Ar B near the center, 
will serve as a viable means of extrapolation beyond the u = 20 surface. 

For each case, we calculated the exponents in Table [5J Here, the energy density for a constant 
r observer will be 

P = ~T\ = — 5- (r 2 v T uu - 2r :U r :V T uv + r 2 u T vv ) , (62) 

where the dominant term in r — > limit will be proportional to 

p«(£) 2 , (63) 

where it is proportional to the electric energy density. Thus, the electric energy density is propor- 
tional to 

p(r) oc ~ r 2 ( B ~ l \ (64) 

Therefore, roughly speaking, if B is less than 1, it will be evidence that the energy density p diverges 
around the center, showing that there is a central singularity. 
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Figure 6: Numerical results for P = while varying from e = 3.5 to e = 4.0. 
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Figure 7: Charge q(r) at v = 180 for each e. 




Figure 8: Ricci scalar for u = 12.5, 15, and 17.5, for the e = 0.2 case. 



e 


A 


B 


0.3 


0.0694 


0.7534 


0.35 


0.0332 


1.0327 


0.38 


0.0203 


1.2035 


0.4 


0.0144 


1.3195 


0.5 


0.0019 


2.0176 



Table 2: Numerical fitting for q(r) to the functional form q(r) = Ar B for each value of e. We fitted 
from u — 19 to u — 20 to determine the behavior around the center. 
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Therefore, we may conclude that, when e = 0.3, there is surely a central singularity while for 
the other cases this is not evident. Of course, via the singularity theorem, if there is an apparent 
horizon, there will be a singularity anyway. For the e = 0.4 and 0.5 cases, however, there is 
no horizon. In these cases, these exponents are circumstantial evidence that there may not be a 
singularity around the center. (This conclusion is consistent with the result of |16j.) Of course, 
if the distribution of the charge drastically changes beyond u = 20, the central structure will be 
modified and our extrapolation will be meaningless. Therefore, our prospect depends, in all respects, 
on circumstantial evidence and not on real evidence. 

4.1.3 Mass inflation 

There are several studies in the literature that closely investigated the mass inflation in three- 
dimensional black holes |21) . One of the interesting predictions was that there is a curvature 
singularity for a certain curvature component, whereas some of the other curvature components do 
not diverge as v increases. 

We observed the Ricci scalar R (Figure [8]) , where 

t~, ^ / T ,UV , QQ,UV & ,V,Q ,V \ //ir-\ 

R = -J 1 ~ , (65 

or \ r a z J 

Tf lu l tl l 1 ' oc T vv /a 2 for an out-going null observer l^ (top in Figure [9]), T^ v n^n v oc T uu /a 2 for an 
in-going null observer n M , for some out- going (middle in Figure [3]) and in- going null slices (bottom 
in Figure© (u = 12.5, 15, and 17.5; v = 120, 150, and 180). Note that we choose ({u,v,ip}) 

I* = —(0,1,0), (66) 
a 

n M = ^(1,0,0), (67) 
a 

so that l^lfj, = n^n^ = and l^n^ = — 1. Thus, any time-like geodesic t M can be decomposed by 

i" = al» + bn^ (68) 

with constraint ab = 1/2. Hence, any time-like observer will see local energy density by 

T^tH v = ^ (b 2 T uu + a 2 T vv ) + ^ . (69) 

Therefore, if there is a signature for a case in which T MU /a 2 or T vv /a 2 increases exponentially, then 
it is the evidence of mass inflation. 

The results arc interesting: there is no evidence that the Ricci scalar or T vv /a 2 diverge, while 
Tuu/cv 2 shows exponential growth. Thus, this provides clear evidence that there is mass inflation and 
that a curvature singularity forms around the inner Cauchy horizon at the v — > oo limit. However, 
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Figure 9: Top and middle: T vv /a 2 and T uu /a 2 for u = 12.5, 15, and 17.5, for the e = 0.2 case. 
Bottom: T uu /a 2 for v = 120, 150, and 180, for the e = 0.2 case. 
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this behavior of curvature functions is very different from that noted in four-dimensional cases. For 
four-dimensional cases, almost all related functions (e.g., Ricci scalar, Kretchmann scalar, Misner- 
Sharp mass function) diverge exponentially, while the Ricci scalar and T vv /a 2 components do not 
diverge in three dimensions. Therefore, this is a typical property of three-dimensional charged black 
holes. 

4.2 Semi-classical back-reactions by Hawking radiation 

In this section, now we turn to the semi-classical effects and see dynamical back-reactions for the 
inside structures. 

4.2.1 Results 

We put P = 0.1 and vary parameters e (Figure [T0|) . as was done in the previous section. 

For the e = 0, 0.075, and e = 0.15 cases, there is a new horizon where r„ = r „ = 0. For the 
e = case, there is a space-like singularity. Hence, between the r u = r „ = horizon and the 
central singularity, there is a fuzzy region and a Cauchy horizon along an ingoing null direction. 
The r. v = and r )U = horizons can be separated as in the e = 0.2 case, while for small charge 
cases, they are mostly foliated. Because r „ and r u both change their signs around the foliated 
region, any observer who moves through the horizon sees an increase of areal radius r. Therefore, 
it is a type of a bottleneck in a one-way wormhole. Note that the size of the bottleneck is on the 
order of < P = 0.1. 

For charged black holes in four dimensions, there was a similar effect. However, in those cases, 
the charge played a crucial role [12], while in three dimensions, the charge is not essential to see 
such a bottleneck, as we can see in the e = case. Rather, it is the effect due to the negative energy 
via semi-classical effects that concentrated around the center. This bottleneck becomes a curvature 
singularity as the curvature functions increase exponentially around these regions. Here, the charge 
will make the bottleneck relatively regular. This phenomenon will be discussed in Section 14.2.21 

As we increase the charge, two horizons are separated (in the e = 0.2 case) and the r u = 
horizon even disappears (in the e = 0.3 case). This suggests that the amount of the negative energy 
decreases as the charge increases, approaching and exceeding the extreme limit (e = 0.4 and 0.5). In 
these cases, the inner apparent horizon can be space-like while the outer horizon remains space-like. 
This implies that the null energy condition is violated only for inside and that it is not violated 
around the outer horizon. This will be discussed in Section 14.2.31 
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A=0.2, A/8tt= -0.0001, P=0.1 



e=0 e = 0.075 




Figure 10: Numerical results for P = 0.1 (with Hawking radiation) while varying e. 



2G 





Figure 11: Ricci scalar for P — 0.1 along u = 15 when e = 0.05, 0.075, 0.1, 0.5, and 0.2 (Upper) and 
along v = 180 when e = 0.05, 0.075, 0.1, 0.5, and 0.2 (Lower). 

4.2.2 Formation of weak curvature singularities 

Figure [Til show Ricci scalar functions for P = 0.1 along u = 15 and v = 180 while varying e. Until 
e = 0.15, the r, v = horizon and the r „ = horizon are foliated; hence, the inner horizon becomes 
a bottleneck of a wormhole, after which the Ricci scalar increases exponentially. Unless we assume 
a large number of scalar fields AT so that R < Af, the internal structure becomes a type of curvature 
singularity with a non-zero area. 

Note that this is not due to mass inflation, because mass inflation does not make the Ricci scalar 
diverge, as shown in Section [4.1.31 Also, this effect exists even for the e = case. Therefore, we 
can clarify that this is a new type of curvature singularity. 

As was noted in Section 15751 the bottleneck occurs via the violation of the null energy condition 
and the concentration of the negative energy. If the effect of the re-normalized energy-momentum 
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tensor is dominant (around the order of r < P) compared to the classical fields, then the region is 
no longer semi-classical. Dvali noticed that such a region is not physically well defined and is hence 
physically problematic. It is interesting that curvature functions become trans-Planckian around 
such dangerous regions. Such a curvature singularity is found by the authors for the first time, 
and we will term it a Dvali curvature singularity to distinguish it from the mass inflation curvature 
singularity. 

Both types of curvature singularities are weak, in the sense that the area is non-zero while the 
curvature becomes trans-Planckian |29] . One can compare the order of magnitude of such curvatures 
between the Dvali curvature singularity and the mass inflation curvature singularity. The former 
is dominant compared to the latter in terms of the order of magnitude. Therefore, the following 
conclusions are reasonable: 

• For charged black holes with sufficiently small charge, a space-like singularity partljn changes 
to a Dvali curvature singularity. 

• For large charged black holes, we can see a mass inflation curvature singularity in the large v 
limit. 

• Between the two limits, both types of curvature singularities can exist. 
4.2.3 Energy-momentum tensors and thermodynamic stability 

Figure [T2l shows the typical distribution of the energy-momentum tensor components of our simula- 
tions. The first plot in Figure [T^] is T vv , showing that it is positive during the gravitational collapse 
< 20). 

After the gravitational collapse, the inside structure is governed by the negative energy-momentum 
tensors T uu and T vv . The second plot in Figure Q21 shows that T uu is negative around the outer 
part of the r M = horizon and changes its sign around the inner part of the r u = horizon. This 
is consistent given that 

r,uu\r, u =o = -SnrT uu , (70) 

and hence, if r tU increases along the ingoing null direction at r u — 0, T uu should be negative; if r u 
decreases along the ingoing direction at r u — 0, T uu should be positive. 



2 If the bottleneck approaches too closely to the outer horizon, and a space-like singularity can appear, as we 
see in the e = case. This is due to the fact that the T vv component is positive around the outer horizon via 
thermodynamic stability, as we shown in the next section. 
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Figure 12: Energy-momentum tensor components T vv , T uu , and T vv around the inner horizon, and 
(T^ v ) around the outer horizon, for e = 0.2 and P = 0.1. 

The third plot in Figure [T^] is T vv around the inner horizon. It is negative around the inner 
apparent horizon. Note that around the inner horizon, 

r,vv\r, v =a = SnrT vv . (71) 

If r „ — horizon is inside of the black hole, to increase r tV from negative to positive, the inner 
horizon should be space-like, and it is equivalent to T vv < 0. Therefore, this is a consistent behavior. 

Finally, the fourth plot in Figure [T2l shows the (T*) components. Of course, the outer apparent 
horizon is space-like or time-like if and only if T vv is positive or negative. However, we know 
that T vv = T^ v + (T^ v ) and that competition exists between the two parts. The first classical 
part is always positive, and hence, to see a time-like horizon, the second quantum part should be 
negative. Moreover, it has to be greater than the first part. However, the fourth plot shows that 
the quantum part (T^ v ) is already positive around the outer horizon. Therefore, there is no hope 
of seeing the evaporation of a black hole in three-dimensional cases. This is consistent with the 
known thermodynamic stability of three-dimensional black holes, and this calculation confirms the 
stability in an independent manner [7]. 

4.2.4 Modifications of causal structures due to semi-classical effects 

In this subsection, let us summarize the results of previous sections. 
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Figure 13: Causal structures with Hawking radiation. Type 1, 2, and 3 in Figure [5] should be 
modified to Type T, 2', and 3'. 

In Section 14.1.11 we observed causal structures when there is no Hawking radiation. In Sec- 
tion 14.2.11 and Section 14.2.31 we included Hawking radiation and observed some interesting modifi- 
cations: 

• Space-like singularities are partly changed to a bottleneck of a wormhole (the Dvali curvature 
singularity). 

• Inner horizons at finite v can be space-like. (This depends on the details of the energy- 
momentum tensors of inside the black hole.) 

• A v —¥ oo inner horizon remains for sufficiently large charge cases. 

• The outer apparent horizon is thermodynamically stable. 

According to Section 14.1.31 and Section 14.2.21 we understand that the v — > oo inner horizons and 
the bottleneck of the wormhole become curvature singularities. The former is due to mass inflation 
and the latter is caused by the concentration of negative energy. 

Therefore, for Type 1 (Figure [SJ, a part of the space-like singularity will be changed by a Dvali 
curvature singularity (left of Type 1', Figure I13p. We cannot trust beyond the Dvali curvature 
singularity, since the curvatures become trans-Planckian. However, as was discussed in Section f3.51 
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if there are sufficiently large numbers of massless fields that contribute to Hawking radiation without 
changing P oc Afh, beyond the r „ = r v = horizon can be trustable in the semi-classical sense, 
as the curvatures, for example the Ricci scalar, will decrease in proportion to 1/JV 2 . If this is true, 
our solution implies a branching-off geometry (right of Type 1', Figure [TBI . For Type 2 (Figure [5]), 
again, the space-like singularity should be changed by a Dvali curvature singularity and the inner 
Cauchy horizon should be changed by the mass inflation singularity (Type 2', Figure H3)). For 
Type 3 (Figure [5]) , the semi-classical effect can modify the inner horizon to be space-like and there 
will be a mass inflation singularity in the v — > oo limit (Type 3', Figure [T3")) . Of course, in that 
limit, the out-going Cauchy horizon will be also a type of mass inflation singularity. In this limit, 
we may not see the Dvali curvature singularity. 

5 Conclusion 

In this paper, we investigated the dynamical formation and evolution of three-dimensional charged 
black holes. As a concluding remark, we summarize our new results and point out some interesting 
issues in terms of two points: classification of causal structures for dynamical black holes and cosmic 
censorship in three dimensions. 

5.1 Causal structures 

We study gravitational collapses of charged matter fields and classified causal structures (Figure [5]): 
Type 1 for an almost neutral case, Type 2 for a small charge, Type 3 for a large charge, and Type 4-1 
or Type 4-2 for an excessive charge. For near-extreme cases, we show evidence that Type 4-2 will be 
more probable. When there is an inner Cauchy horizon in v — > oo limit, a local energy density for an 
in-going null observer increases exponentially, and this confirms mass inflation. This will eventually 
create a curvature singularity in a large v or large u limit. However, this was significantly different 
from the four-dimensional case, as in those cases, other types of curvature functions also diverge. 

As we include semi-classical effects via Hawking radiation, we note that the outer apparent 
horizon is not affected and is thermodynamically stable. However, the internal structure will be 
significantly modified. Negative energy is concentrated near the center, and this can create a 
bottleneck in a wormhole. Such a bottleneck is not due to a charge but is caused by negative energy. 
This neck occurs when the semi-classical effects dominate the classical effects. We observed that 
the region around the bottleneck becomes a curvature singularity (the Dvali curvature singularity) . 
Moreover, the inner horizon can be space-like via negative energy. Therefore, we delineate two 
types of curvature singularities in three dimensions: the Dvali curvature singularity and the mass 
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inflation curvature singularity. Therefore, as we summarize these contents, the causal structures 
Type 1, Type 2, and Type 3 should be modified to Type 1', Type 2', and Type 3' (Figure [BJ . 

5.2 Comments on cosmic censorship 

One interesting finding was related to cosmic censorship. There is some controversy about whether 
or not there is a violation of cosmic censorship in anti de Sitter space |36) |37j . In terms of BTZ 
solutions, there appears to be a time-like or null singularity for M < Q 2 InQ 2 or M = Q = cases. 
Therefore, as long as we rely on the BTZ solution, a violation of weak cosmic censorship is not odd. 

However, if we believe that all realistic black holes should be formed by a gravitational collapse, 
our conclusions seem to say that there is no violation of weak cosmic censorship fSection 14. 1.21) . 
When there is no horizon and excessive charge, the tendency of the charge distribution shows that 
the energy density around the center is regular. 

The absence of a naked singularity has also been observed in four-dimensional gravitational 
collapses [H][12]. However, in that case, as we give an excessive charge, a near-extreme black hole 
was formed rather than a disappearance of a black hole. This is an interesting difference between 
three dimensions and four dimensions. 

In terms of strong cosmic censorship, there can be two ways to touch the Cauchy horizon. First, 
there is the Cauchy horizon along the in-going null direction (e.g., the thin red dotted line in Type 3 
in Figure [5]). Second, there is the Cauchy horizon along the out-going null direction (e.g., thick 
black dashed line in Type 3 in Figure [S]) . Given that this black hole does not evaporate, for the 
former Cauchy horizon, an ingoing observer requires infinite advanced time v (although this does 
not imply infinite proper time). Thus, via mass inflation, it will surely be a curvature singularity. 
Therefore, it is likely not possible to touch in the semi-classical regime. Of course, it may be possible 
for an observer to touch the latter Cauchy horizon. However, we know that when there are two 
horizons (a space- like outer horizon and a time- like inner horizon), it is not evident whether there 
is a time-like singularity for an ingoing null observer, as we noted in Table [5] If this is not the case, 
then strong cosmic censorship appears to be safe. 
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Appendix A. The Birkhoff theorem for charged BTZ black 
holes 

In this appendix, we briefly examine the Birkhoff theorem for charged BTZ black holes. There is 
an argument that the only spherical symmetric solution of (2 + l)-dimensional space-time with a 
negative cosmological constant is a BTZ black hole or the self-dual Coussaert-Henneaux space-time 
|38) . This was applied only for neutral cases, and hence it is necessary to present an argument about 
the uniqueness of the solution when the electric charges are involved. In this appendix, we prove 
that any spherical symmetric charged black holes in a three-dimensional anti de Sitter background 
will approach the charged BTZ asymptotically. 

The metric of a spherically symmetric space-time can always be cast in the form 

dg 2 = _ e 2*(T,t)f, t )dt 2 + -^—dr 2 + r 2 d V 2 (72) 

f(r,t) 

involving the two arbitrary functions ij)(r,t) and f(r,t). We put the gauge field by 

A t (r) = Qln- l , (73) 

as the electric field strength must be of the order 1/r asymptotically. We may put the gauge field 
as an arbitrary function for completeness. However, this restriction of the configuration space is 
reasonable because the conservation of the energy-momentum tensor, V^T^ 1 ,, = 0, constrains the 
form of the gauge field as Equation (jT3|) . 

The independent Einstein equations are as follows: 



i , Q 2 , f' 2 
i 



+ ^e~ 2lp + 4— = 0, (tt)-component (74) 

l z r z 2r 



e~ 2v 7 = 0, (tr) -component (75) 

2lp -l h -ip' = 0, (it)- component (76) 



2rf 2 

P r 2 2r r 



where ' and ' represent the derivatives with respect to r and t, respectively. 

From the (rr)-component and the (tt)-component of the Einstein equations, we have ip = ip{t). 
In addition, a time redefinition can absorb ip(t) and it will make f(r,t) be the only function to be 
determined. From the (ir)-component, we have / = f(r). Finally, from the (tt)-component of the 
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Figure 14: The real part of the scalar field as a function of v around the outer apparent horizon. 
Here, we used e = 0.2 and P = 0. The gradients of each lines are approximately —0.29. 

Einstein equations, we obtain 



where C is the integration constant and we can choose by —M. Therefore, this proves our assertion. 
Any spherical symmetric charged black holes in the three-dimensional anti de Sitter background 
will approach the charged BTZ solution. 

To sure this expectation, we test the dynamics of scalar field around the outer apparent horizon. 
For e = 0.2 and P = 0, we observe the behavior of the scalar field (Figure . For convenience, we 
observed the real part of the scalar field as a function of v around the outer apparent horizon (2.78 < 
u < 2.87); the imaginary part is also qualitatively similar. The scalar field falls off approximately 
~ v~ p with p ~ 0.29. Therefore, in the v — > oo limit, the scalar field will decay to zero. Therefore, 
this makes sure that our solutions will emerge to a vacuum solution in the anti de Sitter. Then, 
our previous proof shows that our results surely emerge to BTZ solutions. 

Appendix B. Consistency and convergence tests 

Now we test consistency and convergence of this paper. 

First, to check the consistency, we used one of constraint equations. As an example, we choose 
the constraint function Equation (|33[) 




(77) 



g >v - 2gd + 2rzz + 4P(d, v - d 2 )\ 



(78) 



\g, v \ + \2gd\ + \2rzz\ + 4P(\d, v \ + \d 2 \) 
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to see the constraint. The value should be exactly zero, but due to numerical errors, it may not 
be zero. The normal scale of the constraint function is on the order of 1; hence, if the value is 
sufficiently less than 1, it is certain that our simulations give a sound result. Top in Figure 1151 
shows the constraint function for the P = 0.1 and e = 0.2 case along the u — 5, 10, and 15 slices. 
The result is sufficiently small and the order is less than 1%. Note that when the constraint becomes 
on the order of 1%, the r tVV actually approaches zero and the behavior of the constraint function 
becomes unstable. However, this is not a significant problem for entire causal structures. 

Second, to check the convergence, we first see \r^ 2 ) ~ r (i)l/ r (2) an d M r (4) — r (2)\/ r (4) along the 
u = 5, 10, and 15 slices, where r n means an n x n finer simulation (Middle in Figure [TS"]) . We 
checked for the e = 0.2 and P = case. For each slice, two curves are almost same; this shows 
the second order convergence. Numerical errors are less than 0.000001%. Next, we test the e = 0.2 
and P = 0.1 case (Bottom in Figure 1X5"]) and compare |r( 2 ) — rn\\/r(2) and 2|r(4) — r(2)|/?"(4). As 
we noted, it shows the first order convergence. However, the error is sufficiently small again, and 
hence it is reasonable to trust our numerical calculations. 
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Figure 15: Top: constraint function \g. v -2gd+2rzz+AP(d iV -d 2 )\/(\g tV \ + \2gd\ + \2rzz\+AP(\d iV \ + 
\d 2 \)) for P = 0.1 and e = 0.2 along u = 5, 10, and 15. Middle and Bottom: convergence tests 
l r (2) — r (i)\/ r (2) and 2|r( 4 ) — r^ 2 )\/r^ for P = and P = 0.1 with e = 0.2 along u — 5, 10, and 15, 
where r„ means an n x n finer simulation. 
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